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Abstract : Multifractional Brownian motion (mBm), denoted here by X, is one 
of the paradigmatic examples of a continuous Gaussian process whose pointwise 
Holder exponent depends on the location. Recall that X can be obtained (see 
e.g. [BJR97, AT05]) by replacing the constant Hurst parameter H in the stan- 
dard wavelet series representation of fractional Brownian motion (fBm) by a smooth 
function H(-) depending on the time variable t. Another natural idea (see [BBCI00]) 
which allows to construct a continuous Gaussian process, denoted by Z, whose point- 
wise Holder exponent does not remain constant all along its trajectory, consists in 
substituting H(k/2 J ) to H in each term of index (j, k) of the standard wavelet series 
representation of fBm. The main goal of our article is to show that X and Z only 
differ by a process R which is smoother than them; this means that they are very 
similar from a fractal geometry point of view. 

Keywords: Fractional Brownian motion, wavelet series expansions, multi- 
fractional Brownian motion, Holder regularity. 



1 Introduction and statement of the main results 

Throughout this article we denote by H(-) an arbitrary function defined on the 
real line and with values in an arbitrary fixed compact interval [a, b] C (0, 1). 
We will always assume that on each compact JC C M, H(-) satisfies a uniform 
Holder condition of order > b i.e. there is a constant c\ > (which a priori 
depends on /C) such that for every t\,t2 £ IC one has, 

iHih) - H(t 2 )\ <c 1 |i 1 -t 2 |' 3 ; (1) 

typically H(-) is a Lipschitz function over M. We will also assume that 
a = mi{H(t) : t e M} and b = sup{H(t) : t E M}. Recall that multi- 
fractional Brownian motion (mBm) of functional parameter H(-), which we 



2 



Antoinc Ayache and Pierre R. Bertrand 



denote by X = {X(t) : t e IR}, is the continuous and nowhere differentiable 
Gaussian process obtained by replacing the Hurst parameter in the harmo- 
nizable representation of fractional Brownian motion (fBm) by the function 
H(-). That is, the process X can be represented for each t E Mas the following 
stochastic integral 

x[t) = lmW^ mi (2) 

where dW is "the Fourier transform" of the real-valued white-noise dW in the 
sense that for any function / e L 2 (M) one has a.s. 

[ f(x)dW(x)= [ M)dW(Z). (3) 
J JR J IR 

Observe that (3) implies that (see [C99, ST06]) the following equality holds 
a.s. for every t, to within a deterministic smooth bounded and non-vanishing 
deterministic function, 



Ie 



\H(t) + l/2 



dW(0 = [ {\t- «|*(*)-V2 - | S |*(*)-V2\ dW ( s) . 



Therefore AT is a real- valued process. MBm was introduced independently in 
[PLV95] and [BJR97] and since then there is an increasing interest in the 
study of multifractional processes, we refer for instance to [FaL, S08] for two 
excellent quite recent articles on this topic. The main three features of mBm 
are the following: 

(a) X reduces to a fBm when the function H(-) is constant. 

(b) Unlike to fBm, ax = {ax{t) : t G IR} the pointwise Holder exponent of 
X may depend on the location and can be prescribed via the functional 
parameter H(-); in fact one has (see [PLV95, BJR97, AT05, AJT07]) a.s. 
for each t, 

a x (t) = H{t). (4) 

Recall that ax the pointwise Holder exponent of an arbitrary continuous 
and nowhere differentiable process X, is defined, for each t € M, as 



f \X(t + h)-X(t)\ J 
a x {t) = supja e M+ : limsupJ— i ^ = 0| 



(5) 



(c) At any point t £ M, there is an fBm of Hurst parameter H(t), which is 
tangent to mBm [BJR97, F02, F03] i.e. for each sequence (p n ) of positive 
real numbers converging to 0, one has, 

l im i aw J X(t + p n u) -X(t) u e R \ law{£ > (u) . u g R ^ (6) 

n— >oo I n ( ' 

K Pn ) 

where the convergence holds in distribution for the topology of uniform 
convergence on compact sets. 
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The main goal of our article is to give a natural wavelet construction of a 
continuous and nowhere differentiable Gaussian process Z = {Z(t)} te jp^ which 
has the same features (a), (b) and (c) as mBm X and which differs from it by 
a smoother stochastic process R = {R{t) : t G R} (see Theorem 1). 

In order to be able to construct Z , first we need to introduce some notation. 
In what follows we denote by {2^ 2 ip(2 : > x — k) : (j, k) G /Z 2 } a Lemarie-Meyer 
wavelet basis of L 2 (R) [LM86] and we define & to be the function , for each 
{x, 9) e Rx R, 

'<*»>- L^W&i* < 7) 

By using the fact that ^ is a compactly supported C°° function vanishing on 
a neighborhood of the origin, it follows that \P is a well-defined C°° function 
satisfying for any (/, m, n) G IN 3 with / > 2, the following localization property 
(see [AT05] for a proof), 

c 2 = sup {2+\x\f\{d™d^){x,6)\ <co, (8) 

6e[a,b],xeM 

where d^dgty denotes the function obtained by differentiating the function 
n times with respect to the variable 9 and m times with respect to the 
variable x. For convenience, let us introduce the Gaussian field B = {B(t, 9) : 
(t, 9) e Rx (0, 1)} defined for each (i, 9) e R x (0, 1) as 

Observe that for every fixed 9, the Gaussian process B(-,6) is an fBm of 
Hurst parameter 6 on the real line. Also observe that mBm X satisfies for 
each t G R, 

X(t) = B(t,H(t)). (10) 

e iti _ i 

By expanding for every fixed (t, 6), the kernel function £ i— > ^ 9+] ^ 2 in the 

orthonormal basis of L 2 {R), {2^ /2 (27r) 1/2 e 42 ^ fc «^(-2^0 : U, k) G % 2 } 
and by using the isometry property of the stochastic integral in (9), it follows 
that 

oo oo 

W)= E E 2-'% Vb {#(2''t-M)-#(-M)}, (ii) 

j=—oo k— — oo 

where {sj,k '■ G Z 2 } is a sequence of independent Af(0, 1) Gaussian 

random variables and where the scries is, for every fixed (t, &), convergent in 
L 2 (fi); throughout this article ft denotes the underlying probability space. In 
fact this series is also convergent in a much stronger sense, see part {%) of the 
following remark. 

Remark 1. The field B has already been introduced and studied in [AT05]; we 
recall some of its useful properties: 
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(i) The series in (11) is a.s. uniformly convergent in (t, 6) on each compact 
subset of IR x (0, 1), so B is a continuous Gaussian field. Moreover, com- 
bining (10) and (11), we deduce the following wavelet expansion of mBm, 

oo oo 

X ^= E E 2- m %,k{n*t-k,H(t))-$(-k,H(t))}. (12) 

j= — oo k— — oo 

(ii) The low frequency component of B, namely the field B = {B(t, 9) : 
(t, 9) £ Mx (0, 1)} defined for all (t, 6) £ R x (0, 1) as 

— 1 OO 

= E 2-^ fc {#(2*t-M)-#(-M)}, (13) 

j= — oo k— — oo 

is a C°° Gaussian field. Therefore (1) and (10) imply that the low frequency 
component of the mBm X, namely the Gaussian process X = {X(t)} te jp> 
defined for each t £ IR as 

— 1 oo 

E E 2- mt) ej,k{*(yt-k,H(t))-V(-k,H(t))}, (14) 

j — — oo k— — oo 

satisfies a uniform Holder condition of order [3 on each compact subset K, 
of IR. Thus, in view of (b) and the assumption b < (3, the pointwise Holder 
exponent of X is only determined by its high frequency component, namely 
the continuous Gaussian process X = {X(t)} tf -jft defined for each t £ IR 
as 

+ oo oc 

*(*) = E E z-^s^^yt-^Hw-n-KHit))). (15) 

j—0 k— — oo 

Definition 1. The process Z = {Z(t) : t £ defined for each t e M as 

oo oo 

Z W= E E 2-''*( fc /*>^ fc {#(2*t-^ 

j — — oo k— — oo 

(16) 

In view of (11) it is clear that the process Z reduces to a fBm when the 
function H(-) is constant; this means that the process Z has the same feature 
(a) as mBm. 

Remark 2. Using the same technics as in [AT05] one can show that: 

(i) The series in (16) is a.s. uniformly convergent in t on each compact interval 
of M; therefore Z is a well-defined continuous Gaussian process. 
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(ii) The low frequency component of the process Z, namely the process Z = 
{Z(t) : £ G R} defined for all* G M as 

— 1 OO 

^(*)= E E ^ H{k/23) Hk{n^t-k,H(k/^))-<P(-k,H(k/^))} 7 

j — — oo k— — oo 

(17) 

is a C°° Gaussian process. The pointwise Holder exponent of Z is therefore 
only determined by its high frequency component, namely the continuous 
Gaussian process Z = {Z(t) : t G M} defined for all £ G M, as 

+ OC oo 

^(*)=E E ^~ jH{k/2S) Sj,k{^C2 j t-k,H(k/2 j ))-^(-k,H(k/2 j ))}. 

j=0 fe=-oo 

(18) 

It is worth noticing that if one replaces in (18) the Holder function H(-) by 
a step function then one recovers the step fractional Brownian motion which 
has been studied in [BBCI00, ABLV07] . 

Let us now state our main result. 

Theorem 1. Let R = {R(t) : £ G M} be the process defined for any t £ 1R as 

R(t) = Z(t) - X{t). (19) 

Let K. be a compact interval included in M. Then, if a and b satisfy the fol- 
lowing condition: 

1 - b > (1 - o)(l - ab- 1 ), (20) 

there exists an exponent d G (6, 1] , such that the process R satisfies a uniform 
Holder condition of order d on fC. More precisely, there is f2* an event of 
probability 1, such that, for all u G f2* and for each (£o,£i) G K? , one has 

\R(ti,u) - R(to,u)\ < C 1 (uj)\t 1 -t \ d , (21) 

where C\ is a nonnegative random variable of finite moment of every order 
only depending on Q* and K. 

Remark 3. We do not know whether Theorem 1 remains valid when Condition 
(20) does not hold. Figure 1 below indicates the region V in the unit cube 
satisfying (20). 

Thanks to the previous theorem we can obtain the following result which 
shows that Z and X are very similar from a fractal geometry point of view. 

Corollary 1. Assume that a and b satisfy (20), then the process Z has the 
same features (a), (b) and (c) as mBm. 

Throughout this article, we use [x] to denote the integer part of a real 
number x. Positive deterministic constants will be numbered as c\, C2, . . . while 
positive random constants will be numbered as C\,C<2,, . . .. 
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Fig. 1. the region T> in the unit cube satisfying (20) 



2 The main ideas of the proofs 

In the reminder of our article we always assume that Condition (20) is satist- 
fied and that diam(/C) := sup{|u — v\ : (it, v) 6 JC] < 1/4. Also notice that we 
will frequently make use of the inequality 

log(3 + x + y) < log(3 + x) x log(3 + y) for aU (a;, y) 6 R\. (22) 

Let us now present the main ideas behind the proof of Theorem 1. Firstly 
we need to state the following lemma which allows to conveniently bound the 
random variables Ejj,- It is a classical result we refer for example to [MST99] 
or [AT03] for its proof. 

Lemma 1. [MST99, AT03J There are an event Q* of probability 1 and a 
nonnegative random variable C2 of finite moment of every order such the 
inequality 

< C 2 MVlog(3+|j| + |A;|), (23) 
holds for all u> £ fl* and j, k G Z. 

Proof of Theorem 1. In view of Remark 1 (ii) and of Remark 2 it is sufficient 
to prove that Theorem 1 holds when the process R is replaced by its high 
frequency component, namely the process R = {R{t) : t G ]R} defined for 
each t G 1R as 

R(t) = Z(t) - X(t). (24) 
Let g^k be the function defined on 1R x M, by 



g jtk (t, 6) = 2-^i^{2H -k,6)- !?(-*, 9)}. (25) 
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It follows from (24), (15), (18), (25) and (23) that for any uj £ Q* , 

\R{t u u)-R{t ,u)\<C 2 {w)Y, E V^g(3 + j + \k\) 

j—0 k— — oo 



(26) 



9j> fc (ti, H(k/V)) - g jtk (t , H (fc/2*)) - g jt k {t u H(h)) + g jt k (t Q ,H(t )) 

Next, we expand the term <7j,fc(tj, H(t)) with i — or 1 and r = fi or fc/2 J 
with respect to the second variable in the neighborhood of H(to). Indeed, 
since the function & is C°°, the functions gj. k are also C°°; thus we can use 
Taylor-Lagrange Formula of order 1 with an integral reminder and we get 

g^kituH^)) = g jtk {h,H{t )) + (H(h) - H(t )){de 9j . k ){tuH{t )) (27) 

+ (H(h) - if (t )) 2 / (1 - T^gg^tuHito) + T(H(h) - H(t ))) dr, 
Jo 

g jtk (t ,H(k/2i)) = g jik (t ,H(t )) + {H(k/V) - H(t ))(d g g 3 , k )(t , H(t )) (28) 

+ (H(k/V) - H(t )) 2 f (1 - T)(d 2 e9jtk )(t , H(t ) + T(H(k/V) - H(t ))) dr, 
Jo 

and 

g jtk (t u H(k/2^) = g^tuHito)) + (H(k/2^) - H(t ))(d egj , k )(h (29) 

+ (H(k/2r)-H(to)) 2 f (l-T)(d 2 g g 3 , k )(t u H(t ) + T(H(k/y)-H(t )))dT. 
Jo 

By adding or subtracting relations (27), (28) and (29) the constant terms 
disappear and we get the following upper bound 

9j , k (h,H(k/2i)) -g j<k (t ,H(k/V)) - 9j , k (t u H{h)) + g jik {t ,H{t ))\ (30) 
< \H{H) - H{ta)\ | {degjMu H (t )) 

pi 

+ \H(t 1 )-H(t )\ 2 / (1-r) (d 2 egj , k )(t u H(t ) +r(H(t 1 ) - H(t ))) dr 
Jo 

+\H(k/2?) -H(t )\\(d e g j<k )(t u H(t )) - (d e g j<k )(t ,H(t )) 

+ \H(k/2>) - H(t )\ 2 [ (1 - t) (d 2 e g j!k )(h,H(t )+T(H(k/2i) - H(t ))) 
Jo 



-($&,*)(*(>, #(*<>) + T(H{k/V) - (t ))) 



dr. 



Then, we substitute the previous bound (30) into the inequality (26). We stress 
that the quantities \H(ti) — H(t )\ and l-ff(ti) — H(t n )\ can be factorized 
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outside the sum whereas the quantities \H(k/2 : >) — H(t )\ and \H{k/V) — 
H(to)\ 2 remain inside the sum. We obtain 

\R(h,Lu) - R(to,f>)\ < C 2 (uj)\H(t 1 ) - H(t Q )\ 

+ OO +OO 

x{]T V^g(3 + j + \k\) x \(dg gj , k )(h,H(t ))\} 

j—0 k— — oc 

l-oo +00 



+ C 2 (u)\H(t 1 )-H{t )\ 2 x Yl V^g(3 + j + \k\) 

j—0 k— — oo 

x J (1 - r)|(^ 5j -, fc )(ti,ff(to) + T(H(h) - H(to)))\ dr} 

+ OO +OO 

+ C 2 (cv)x{J2 E A/log(3 + j + |fc|)|ff(fc/2^)-^(*o)l 

j—0 fc— — 00 

x|(9^,fc)(ti,H(to)) - (d e g 3 ,k)(t ,H{t a ))\} 
+ C 2 Hx{^ ]T y/log(3 + j + \k\)\H(k/V)-H(t )\ 2 



j—0 k— — oo 
"1 

}2 



x / (1 - t) Ud 2 e9j , k )(h,H(t ) + r(H(k/V)-H(t ))) 
Jo 1 

-(d 2 e9jtk )(t ,H(to)+T(H(k/2i)-H(t )))\dT}. 

Then using the following two lemmas whose proofs will be given soon, we get 
that 

\R(h,w) -R{to,w)\ < C 2 (u){\H(t 1 )-H(to)\Ai(lC;a,b) + ... (31) 
+ \H(h) - H(t Q )\ 2 A 2 (IC;a,b) + . . . 



+ 



\h -t \ dl Qi{1C;a,b,di) + \t 1 - t \ d2 g 2 ()C;a,b,d 2 )y 



Finally, in view of (1) the latter inequality implies that Theorem 1 holds. M 
Lemma 2. For every integer n > and (t, 9) G JR x (0, +00) one sets 
+00 +00 

A n (t,6):=J2 E mg 3 , k (t,0)W^g(3+3 + \k\). (32) 

j—0 k— — co 

Then one has 

An{fC;a,b) := sup ^A n (t, 9) : (t, 9) G JC x [a, 6] | < 00. (33) 
Lemma 3. For every integer n > 1 and (to,ti,6) G JR 2 x (0, +00) one sets 
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+oc +00 

G n (t o ,t u 0):=J2 E \ H W 2j ) - ff (*o)l" x Vlog(3 +j + \k\) 

j—0 k— — oo 

x|(^,fe)(2 J 'ii -k,8)- (dZg 3 , k )(2Ho -k,6). 
Then, for every integer n>l, there is an exponent d n E (b, 1] such that 

g n (JC;a,d n ) := sup \h - t Q \- d "G n (t Q , t u 6) < oo. (34) 

(t ,t 1 M)elC 2 x[a,b] 

Proof of Lemma 2. From Lemma 4 given in next section, one can deduce 

n +oc +oc 

A„M)<£C£|log2p>£ E J p 2- J 'Vlog(3 + j + |fe|) (35) 

x { I {d n e - p *) {2H-k,0)\ + \ {d n e - p *)(-k, 0)|}. 

Note that the deepest bracket { | (<9™~ p !^) (2^ - M)| + |(5?~ p ^)(-fe,fl)|} 
contains two terms: the first | (dg ~ Px P){2H — k,9)\ depends on t E JC whilst 
the second \{dg~ Pl ^){— k,9)\ no longer depends on t. Therefore, it suffices 
to obtain a bound of the supremum for t E K, of the sum corresponding 
to the first term, then to use it in the special case /C = {0} to bound the 
sum corresponding to the second term. Let us remark that there exists a 
real K > such that K, C [— K, K}. Thus, without any restriction, we can 
suppose that K. = [—K,K]. Next, using (8), the convention that 0° = 1, the 
change of variable k = k' + [2H], the fact that \t\ < K, (22) and the fact that 
z = 2H — [2H] E [0, 1] , one has the following estimates for each p E {0, . . . , n} 
and (t,9) E [-K, K] x [a, b]: 

E E f 2 ^ 8 + J + 1*1) \{d n e - p *){2H -k,6)\ 

j—0 k— — oo 

+ oo +00 

<^EE J p 2^Vlog(3 + j + |fc|)-(2+|2^-fc|)- £ 

j—0 k— — oo 

+ oo +oo 

<c 2 E E J P 2^'Vlog(3 + j + W\ + ¥K) ■ (2 + \2H - [2H] - k'\)~ e 

j—0 k' — — oo 

+ oo 

< c 2 c 3 E f 2 ~ j Vlog(3 + j + 2i K) < oo, (36) 

1=0 

where 

c 3 = sup { (2+|z-A:|)-Vlog(3+|fe|) : ze[0,l]}<oo. (37) 

k— — oo 
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Clearly, (36) combined with (35) implies that (33) holds. ■ 

Proof of Lemma 3. The proof is very technical so let us first explain the 
main ideas behind it. For the sake of simplicity, we make the change of no- 
tation t\ = t + h. Then we split the set of indices {(j, k) 6 IN x Z} into 
three disjoint subsets: V a neighborhood of radius r about t , a subset W 
corresponding to the the low frequency (j < j\) outside the neighborhood V 
and a subset W c corresponding to the the high frequency (J > ji) outside 
the neighborhood V (the "good" choices of the radius r and of the cutting 
frequency j\ will be clarified soon). Thus the sum through which G n (to, t%, 0) 
is defined (see the statement of Lemma 3) can be decomposed into three parts: 
a sum over V, a sum over W and a sum over W c ; they respectively be denoted 
-Bi,n(io, h, 9), B2, n (to, h, 9) and Bz, n {to, h, 9). In order to be able to show that, 
to within a constant, each of these three quantities is upper bounded by \h\ dn 
for some exponent d n > b, we need to conveniently choose the radius r of the 
neighborhood V as well as the cutting frequency ji . The most natural choice 
is to take r = \h\ and 2~ J1 ~ \h\. However a careful inspection of the proof 
of Lemma 7 shows this does not work basically because 2^\h\ does not go 
to infinity when \h\ tends 0. Roughly speaking, to overcome this difficulty we 
have taken r = \h\ v and 2~ J1 ~ |/i| 7 where < 77 < 7 < 1 are two parameters 
(the "good" choices of these parameters will be clarified soon) as shown by 
the following Figure 





t a -\h\i t to + fi t a +\h\" t 



Fig. 2, the three Lemmas and the corresponding subset of indices 
More precisely, j\ is the unique nonnegative integer satisfying 

2- jW < < 2~ J1 (38) 
and the sets V, W and W c are defined by: 

V(t ,h,r))={(j,k)£lNxZ: \k/V - t \ < H"}, (39) 
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V^^jHI^lelxI: \k/V-to\>\hy>}, (40) 
W(t ,h,r 1 , 1 ) = {(j,k) e V c (t ,K V ) : < j < jr} (41) 



and 



W c (i ,/i,r?,7) = {0', fc) e V c (t ,h, V ) : j > j t + 1}. (42) 
It follows from Lemmas 4 to 7 that 

3 

G n (to,to + h,8) = B mtn (to,h, 0) 

m— 1 

n 3 

< E E CZ(log2rB m (t ,h,d) (43) 

p=0 m— 1 



< 



C A (\h\ a+r 'P + |^|(l-7)+7a + | /l |( 7 -r,)(^-l-e)+7^ l g T1+ 1 /2 ( 1 / 1 ^ | ) , 



where the constant C4 = max{c 5 , eg, cio} X)p=o Cn(k ) g2) p does not depend on 
(to, h, 9). In view of (43) and the inequality > b as well as the fact that e is 
arbitrarily small, for proving that (34) holds its sufficient to show that there 
exist two reals < rj < 7 < 1 and an integer I > 2 satisfying the following 
inequalities 

a + Vjb > b 
(1 - 7) + 7a > b 
(7-77)^- 1) + 7 a > b. 

This is clearly the case. In fact, thanks to (20) we can easily show that the 
first two inequalities have common solutions; moreover each of their common 
solutions is also a solution of the third inequality provided that £ is big enough. 
■ 

Before ending this section let us prove that Corollary 1 holds. 

Proof of Corollary 1. Let us first show that Z has the same feature (b) as 
mBm. In view of Theorem 1 and Remark 3 it is clear that the pointwisc 
Holder exponent of R, satisfies a.s. for all t £ M, 

a R (t) > d. (44) 

Next putting together (44), the fact that d > b, (19) and (4) it follows that 
a.s. for all t e 1, 

a z (t) = H(t). 

Let us now show that Z has the same feature (c) as mBm. Let (p n ) be an 
arbitrary sequence of positive real numbers converging to 0. In view of (19) 
and (6), to prove that for each t € M one has 

Um law J Z ( t + Pr^) -Z(t) . u g R 1 = law{B (m) . ueJR}j (45) 

( Pn ) 
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in the sense of finite dimensional distribution, it is sufficient to prove that for 
any u £ IR one has 



Inn E \(W + PnU)-W)Y 



V H(t) ) 

Pn 



0. (46) 



Observe that for all n big enough one has p n \u\ < 1. Therefore, taking K, = 
[t — 1, t + 1] in Theorem 1, it follows that for n big enough, 

E{ c lt+ ji; m f}^ d - Ht,n ^) (47) 

and the latter inequality clearly implies that (46) holds. To have in (45) the 
convergence in distribution for the topology of the uniform convergence on 
compact sets it is sufficient to show that for any positive real L, the sequence 
of continuous Gaussian processes, 

is tight. This tightness result can be obtained (see [B68]) by proving that there 
exists a constant cn > only depending on L and t such that for all n £ IN 
and each m, u 2 £ [-L, L] one has 



if 



j f Z(t + p nUl ) -Z(t) Z(t + p n u 2 ) -Z(t) ^\ 2H(t) 

E <\ -H® -Hit) J < C 1V|«1-«2| (l - 

Pn Pn I 



(48) 

Without loss of generality we may assume that for every n £ IN, p n £ (0, 1]. 
Then by using the fact that (48) is satisfied when Z is replaced by X (see 
[BCI98] Proposition 2) as well as the fact that it is also satisfied when Z is 
replaced by R (this can be done similarly to (47)), one can establish that this 
inequality holds. ■ 



3 Some technical Lemmas 

Lemma 4. For every integer n > and any (t,6) £ IR x IR one has 

(de9j,k)(t,0) (49) 
= j^CH-j log 2Y2-»{{8$-*)9{Vt- M) - 

Proof of Lemma 4- The lemma can easily be obtained by applying the Leibniz 
formula for the nth derivative of a product of two functions. ■ 
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Lemma 5. For each integer n > 1 and (i , h, 9) e M x Ft x (0, +oo) set 
B hn , p (t ,h,6) := Yl j p 2- j9 \H(t ) - H(k/2 j )\ n x 0og(3 + j + \k\) 

(j,k)eV(t ,h,ri) 

x I (d r g l - p V) (2* (t + h)-k,6)- (dl l - p V) (2H -k,6) 

where V(to, h, rj) is the set defined by (39). Then, for all real K > and every 
integers n>l and < p < n, one has 

c 5 := sup \h\- a - nr ' f3 log- p - 1/2 (l/\h\)B hn , p (t , h, d) < oo. 

(to,e)e[-K,K]x[a,b],\h\<l/4 

(50) 

Proof of Lemma 5. It follows from (1) and (39) that 

]T f2^ e \H(t ) - H{k/V)\ n x Vlog(3 + j + |A:|) (51) 

(j,k)eV(t ,h,r,) 

x |(^- p !?)(2»'(to + ft) - k,6) - {d;- p ^)(2H - k, , 
U,k)ev(t Ji,v) 

x\(d^- p ^)(2 j (t +h)-k,0)- (d%- p *)(2 j t -k,6). 

Now let jo > 2 be the unique integer such that 

2~ j °- 1 < \h\ < 2-3°. (52) 

By using (8), (52), the change of variable k = k' + [2 J y], (22) and the fact 
that |to | < K, we can deduce that for any y G [to — 1, to + 1], 

+oo + OC 

^ X! ?2- ja m- p v)V j y-k,e)Wiog(3 + j + \k\) 

j=j +l &:=-oo 

+ oo +00 

^ C2 E E i P2 " j0 ( 2 + l 2 ^ - [ 2j y] -fc'|)-Vlog(3+J+ 2i(|i„|+l)+|fc'|) 

j=jo+l fc' = — oo 

+ OC 

<c 2 -C3 X! j p 2- J 'Vlog(3 + j + 2J(^ + l)) 

<c 6 H a log p+1/2 (l/H), (53) 

where C3 is the constant defined by (37) and the last inequality (in which ce 
is a constant non depending on (to, h, 77)) follows from (52) and some classical 
and easy calculations. 

On the other hand, by using the Mean-value Theorem applied to the 
function dg~ p gj^ with respect to the first variable, (8), the fact that for all 
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2i\h\ < 1 for all j e {0, . . .,j }, (52), (22), (37) and the inequality |t | < K, 
we get that 

jo +oc 

E E i P2 ~ ja x v/ioi(3T7TM 

j— fc— — oc 

\a;- p ^)(v(to + h)-k,e)- (d%- p *)(2h - k, e) 



Jo +oo 



<<*WE E J p 2 j(1 - o) (l + |2 J '*o - [2 J 'to] - fc|)-Vlog(4 + .7 + 2i\t a \ + \k\) 

j—0 k— — oo 

Jo 

< c 2 c 3 \h\ E/2^- a Vlog(4 + .? + 2i K) 

3=0 

<c 7 \h\ a log p+1 / 2 (l/\h\), (54) 

where the constant C7 does not depend on (to,h,rj). Finally, by combining 
(51) with (53) and (54), one can deduce (50). ■ 

Lemma 6. For any (t, h,6) £ 1R x M x (0, +00) and for any integers n > 1 
and < p < n set 

B 2 ,n, P (to,h,e) := E i p 2-^|tf(t ) - H(k/V)\ n x Vlog(3 + j + [fcj) 

C7,fc)eW(t ,/i,»7,7) 

x \(d r e l - p V)(23(t + h)-k,6)- (d%- p V)(2H a - k, 9) 

where W(to, h 7 n,j) is the set defined by (41)- Then, for any real K > 0, one 
has that 

c 8 = sup |/r (1 - 7) - 7 °iog- p - l/2 (i/N)iw*o,M) 

(t ,e)e[-K,K]x[a,b],\h\<l/4 

< 00. (55) 

Proof of Lemma 6. To begin with, note that for any pair of real numbers 
(0 o ,6»i) e (0, l) 2 , one has \6i - 9 \ < 1. Therefore, 

for all (j, k) G JV x ^, |fl-(t ) - H(k/2 j )\ n < 1. (56) 

By using the Mean- value Theorem applied to the function dg~ p \P with respect 
to the first variable combined with (8), we get for all to £ /C and h £ 1R 

(d%- p *)(2i(t + h)-k,6)- (d%- p *)(2H - k, 6) 

< c 2 2 3 \h\ (2 + \2H -k + 2 j uh\)~ e 

for a real number u e (0,1). On the other hand it follows the inequality 
2 J |/i| < 1 for all j e {0, . . . , j\} (which is a consequence of (38)) and from 
triangle inequality that \2H — k + 2^uh\ > \2H — k\ — 1. Therefore 



A process very similar to multifractional Brownian motion 



15 



(d2- p V)(V(t + h)-k,6)- (d^)(2H - k, i 

<c 2 V\h\(l + \2H Q -k\y l 
and as a consequence, we obtain for all (t,h,6) eiRxffix (0, 1) 

jl +oc 

B 2 ,n, P (to,h,e) <c 2 \h\J2 E i P2j(1 " e) (l + fih - fc|)-Vlog(3 + j + 

j— /c— — oo 

Next, making the change of variable k — k 1 + [2Hq] and using triangle in- 
equality as well as the inequality 9 > a, we deduce that for all (t, h, 9) G 
[-#,#] x [-1/4,1/4] x [a, b] 

jl +oc 

5 2l n, P (*o,M) <c 2 |/i|E E ^-^(l + ^o-^'iol-fc'l)-' 

J=0 k' — — oo 

x v /log(3 + j+2i|i | + |A:'|) 



+ 00 



< 



c 2 \h\\ E Vlog(3+|fc'|)(l + |2^ -[2^ ]-fc'| 



» k' — — oo 

E i P 2^' (1 -°) v/log(4 + . 7 +2^|) 

where the last inequality follows from |to| < K and the inequality (22). Set 
z = 2H — [2Ho], obviously z E [0, 1), thus (37) and the latter inequality imply 
that 

B 2 , n , p (t ,h,9)<c 2 -c 3 \h\ \ E.?' P2 ' (1_a Vlog(4 + . 7 + 2> K\) \ , 

J=o 

for all (t, h, 9) e [-K, K] x [-1/4, 1/4] x [a, b]. Finally, in view of the inequal- 
ities 2- 71 < and j\ < log(l/|/i|) (these inequalities are a consequences of 
(38)), we get 

iW*o,M) < c 9 |/.|( 1 -^ a log p+1 / 2 (l/|/ l |), (57) 

where the constant eg does not depend on (to ,h,9). This finishes the proof of 
Lemma 6. ■ 

Lemma 7. For any (t,h,0) £ M x M x [a, 6] and any integers n > 1 and 
< p < n set 

B 3 ,n, P (to,h,e) := E j p 2- J ' fl |ff(t ) -ff(fc/2 J ')| n x Vlog(3 + j + |A:|) 

(j,fe)eW c (t ,/i,77,7) 

x I ~ P ^) ( 2 ' (*o + h) - k,0) - (dg- p V) (2H -k,9) 
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where W c (t , h,rj,j) is the set defined by (42). Then, for every real K > 0, 
for each arbitrarily small real e > and all integer I > 2, one has c w < oo 
where 

c 10 := sup |ft|- (7 -" )(, - 1 - e) - 70 log- p - 1 / 2 (|ft|- 1 )B2,n, P (to,fe,(?) 

(t ,6)e[-K,K]x[a,b],\h\<l/A 

Proof of Lemma 7. By using the triangle inequality combined with (42) and 
(40), one gets, for all (j, k) E W c (t , h, jj, 7), 

I (to + h)- k2~ j \ > |*o - k2~ j + h\> \h\ n - \h\ > cnl^l", (58) 

where the constant en = 1 — 4 r '~ 1 . This means that the integer k necessarily 
satisfies 

\2 j (t + h)-k\ >cii2 J >|". (59) 

In view of (59), let us consider and Tr the sets of positive real numbers 
defined by 



3 

and 



T+ = {\2 3 \t +h)-k\: k G ^ and 2^t + h) - k > Cll 2i\h\^} 

Tf = {\2 j {t + h)-k\: k e Z and k - 2^(t + h) > Cll 2^\h\^)}. 

y fixed j, the set 1 ] 
^ T 3 + ,q) q eIN satisfying for all qE IN, 



For every fixed j, the set T + can be viewed as a strictly increasing sequence 



q + Cli y\h\ r > < T+ q <q + l + Cli y\h\\ (60) 



Similarly, for every fixed j, the set 7^ can be viewed as a strictly increasing 
sequence ( T j~ q ) qe ]N satisfying for all q E IN, 

q + cnVW < rj q < q + 1 + cn2»>|". (61) 

Next, setting 7J- = 7) + U 7^~, it follows that from (8), the triangle inequality, 
the inequality \t + h\ < K + 1, (42), (59), (60) and (61) that 

]T f2-i e \(d2- p V)(2i(t Q + h) - k,e)\y/log(3+j + \k\) 

{].k)eW<=(toM, v ,f) 

<c 2 f2-i e (2+\y(t a + h)-k\y e x ... 

x Vlog(3 + j + \23{t + h)-k\ + 23(K + l)) 

+°° . . -e 

<c 2 ^ ^j p 2-^(2+|r|J y/log(3 + j + \T\+2i(K + l)) 
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+00 +00 



<2c 2 ^ ^jf2-^(2 + (? +c 11 2 J >r) x... 



j=ji+i 9=0 



x A /log(4 + j + q + 2i\h\v + 2i {K + 1)) 



Then, one can use the inequality (2 + x)~ e y / log(4 + x) < ci 2 (2 + x)~ £+e which 
is valid for all nonnegative real number a; where £ is a fixed arbitrarily small 
positive real number and C12 is a constant only depending on e. By combining 
this inequality with (22), (38) and < 1/4, we get 

2 ?T» J (a e "- p f ) (2 j (to +h)-k,6)\ v/log(3 

(j,fc)eW c (to,/i,r/,7) 
+00 +00 

^ 2c 2 E E^ 2 ^'Vlog(3 + .7 + + 1)) x . . . 

j'=j'i+i 9=0 

x (2+ q+ c n 2 j \h\^ _ Vlog(4 + g + cii2J>|") 

<2c 2 c 12 £ 2"^Vlog(3 +j + 20 (K + 1)) x . . . 
j=ji+i 

X (i + ( 1+ y +Cll2j l /l l' 7 )" £+£ ^) 

<c 13 ^ f 2 ~ j Vlog(3 + j + 2l (X + 1))) (l + c n 2 j |/i| r ') ' 
i=ii+i 

+00 

< c^l/ir^" 1 -^ 51 j p 2-^ a+£ - 1 - £ Vlog(3 + j + 2J'(ir + 1)) 

j=ji+i 

< ci5 |/i| (7 -" )( ^ 1 - £)+7a log p+1/2 (l/|^|), 

where ci 2 , C13, C14 and C15 are constants which do not depend on (t n ,h,9). 
Finally, using the latter inequality as well as the triangle inequality and the 
fact that H(-) is with values in [a, b] we get the lemma. ■ 

Acknowledgments. The authors thank the anonymous referee for many 
useful remarks which greatly helped them to improve the earlier version of 
this article. 
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